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FOREWARD 


This  paper  was  presented  at  the  AIAA  1 8th  Aerospace  Sciences 
Meeting,  held  at  Pasadena,  California  on  January  1 4-1 6 ,  1980.  The 
text  has  appeared  as  a  preprint  A1AA-80-0342.  Since  reprints  of  the 
work  are  unavailable,  the  paper,  with  up-dated  errata  sheets  included, 
is  distributed  here  as  a  University  of  Southern  California  School  of 
Engineering,  Department  of  Aerospace  Engineering  Report,  USCAE  138. 

The  work  presents  a  comparison  of  analyses  based  on  the  asymptotic 
theory  for  high  aspect  ratio  wings  with  corresponding  results  from  a 
3_D  ful I -potent ial  computer  code  (FLO  22)  for  oblique  wings  as  well  as 
a  more  conventional  swept  wing  involving 
component  flows. 

The  paper  also  gives  a  survey  and  critique  of  the  relatively  recent 
developments  in  the  lifting-line  theory.  The  material  presented  obvi¬ 
ously  cannot  be  accomodated  by  the  length  of  a  single  standard  journal 
article;  various  parts  of  this  work  are  not  expected  to  reach  journal 
publication  stage  for  some  time.  The  distribution  of  this  report  should 
therefore  serve  a  useful  purpose. 
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ABSTRACT 

Asymptotic  theories  for  high-aspect-ratio  _ 

wings  in  transonic  flow  developed  recently  for 
straight  unyawed  wings  and  for  oblique  wings  show 
that  the  three-dimensional  (3*0)  mixed-flow  calcu¬ 
lations  may  be  reduced  to  solving  a  set  of  2-D 
problems  at  each  span  station.  For  wings  with 
surfaces  generated  from  a  single  airfoil  shape, 
local  slmilutude  exists  In  the  3-0  flow  structure, 
permitting  the  problems  to  be  solved  once  for  al 1 
span  stations.  This  paper  reviews  this  thaoret leal 
development  and  the  related  computational  studies. 

The  essential  elements  in  the  theory  will  be  I  dent r 
fled,  their  roles  are  explained;  their  relationship 
to  the  lifting-line  theory  and  related  classical 
iMthods  are  discussed.  Differences  among  the 
pivoted  (oblique)  wing,  the  swept-back  wing  and  the 
swept-forward-wing  in  the  induced  upwash  are  brought 
out.  Examples  of  similarity  solutions  are  demon¬ 
strated  for  high  subcritical  and  slightly  super¬ 
critical  component  flows;  comparisons  with 
relaxation  solutions  to  a  full  potential  equation 
are  made.  The  examples  include  both  oblique  and 
symmetric  swept  wings,  the  study  also  examine  the 
adequacy  of  the  existing  ful 1 -potential  computer 
code.  Outstanding  problems  remaining  for  subse¬ 
quent  development  are  discussed. 

1.  INTRODUCTION 

Our  sound  understanding  in  the  aerodynamics  of 
wing_sweep,_and  its  use  to  control  the  compress Ibl t*  - 
Ity  effecj,  has  been  limited  in  the  past  mostly  'to  . 
problems  in  the  linear  flow  regimes^ ' "5)  Recent 
interest  in  aircraft  wing  design  has  focused  on  the 
potentiality  for  adopting,  or  modifying,  two- 
dimensional  (2-D)  supercritical  airfoil  data  for 
3*0  eppl i cat  ions. (6-1 I)  The  flow  problem  In  this 
domain  is  necessarily  nonlinear  with  a  transonic 
component  flow.  In  this  paper,  we  review  theore¬ 
tical  and  numerical  studies  on  high-aspect-ratio 
swept  wings  for  this  type  of  flow,  including  sosm 
.very  recent  warkA. - - - - - 
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In  the  domain  of  interest,  the  flow  'laid  far 
from  the  wing  section  should  pertain  to  a  high  sub¬ 
sonic,  or  linear  sonic,  outer  flow,  which  is 
representable. in  the  leading  approximation  by  a 
solution  to  the  Rrandtl-Glauert  equation  correspono- 
ing  to  a  swept  lifting  line  (cf.  sketches  In  Fig.  1',. 
Thus,  the  theoretical,  treatments'***'?)  mentioned 
must  find  roots  in  Prandtl's  original  lifting-line 
ideal’?),  even  though  the  corresponding  inner  pro¬ 
blem  to  be  solved  is  basically  nonlinear  and  may 
Involve  an  imbedded  supersonic  reaion.  Apart  f-om 
gaining  a  better  physical  insight  and  a  greater 
simplicity  in  the  swept-wing  analysis,  the  approecn 
represented  by  Refs.  12-18  reflects  a  desire  for 
implementing  the  current  computer-oriented  3*0 
methods (20-26)whlch,  though  oeinq  very  powerful, 
share  the  common  problem  of  trading  accuracy  witn 
computer  cost  and  storage  requirements,  as  is  well 
known.  While  the  prospect  of  the  next-generation 
computers  with  an  equally  powerful  code  development 
Is  well  In  sight’2®),  they  cannot  directly  answer 
the  question  in  the  aerodynamic  theory:  namely.  In 
what  manner  the  2-D  component  flow  on  a  sweet  wing 
is  approached,  as  the  wing  aspect  ratio  increases. 
The  least  contribution,  to  wnich  the  approach  may 
lead,  is  providing  an  asymptotic  basis  on  wnich  tne 
adequacy  of  the  current  3-D  codes  (their  algorithm 
and  some  of  thair  implicit  assumptions)  can  be  ex¬ 
amined. 

Detailed  expositions  of_the  asymptotic  theorv 
underlying  this  develoomentcan  be  found  in  Refs. 
13-18;  the  following  will  address  more  fully  on 
aspects  of  interest  to  the  aerodynamic  understand¬ 
ing, and  on  examining  its  potentiality  towards 
developing  less  costly  computation  procedures  for 
aerodynamic  design  studies. 

The  following  sections  will  examine  the  esser- 
tlal  elements  of  the  theory,  their  roles  in  control¬ 
ling  the  induced  upwash  field  and  other  solution 
properties.  The  discussion  will  bring  out  a  numoer 
of  aerodynamic. features  distinguishing  oblique 
(pivoted)  wings,  the  swoot-forwara  winqs.  conven¬ 
tions!  swept  winqs;  the  delineation  wl'l  maxe  the 
distinctions  from  classical  lifting-line  anelyses 
apparent.  An  exposition  of  the  reduced  problem  0f 
determining  the  3-D  correction  for  the  (inner)  tran¬ 
sonic  component  flow  is  given  in  section  3  wnere 
some  of  the  subtle  differences  among  Refs.  13-13 
will  be  noted,  and  the  local  similarity  in  the  3-C 
flow  structure  at  well  as  an  unsteady  analogy  are 
explained.  The  features  of  a  recent  development 
based  on  the  ful I -potential  equations  wi'l  also  be 
brought  out.  As  in  most  asymptotic  analysis  o*  this 
type,  the  theory  breaks  down  local lv  in  a  numoer  of 
regions,  of  which  the  most  importent  it  oerheos  the 
vicinity  of  the  apex  of  a  swept  wing  where  the 


canter  tin*  hat  a  slope  discontinuity;  that#  and 
othar  limitations,  as  wall  as  tha  fundamantai 
quad t ion  of  perturbing  a  supercritical  shock- fraa 
solution  ara  examined  in  Section  4.  Examelas  of 
solutions  arc  demonstrated  and  eomparad  with  tha 
ful l-potantia)  results  In  Section  5,  where  soma 
limitation  of  tha  existing  3*0  full-potential  com¬ 
puter  codas  and  a  clarification  of  their  formula¬ 
tion  are  also  given. 


2.  THE  ESSENTIAL  FEATURES  OF  THE  THEORY 


The  basic  idea  used,  in  common  with  Frandtl's 
lifting-line  theory ('“)  for  high-aspect-ratio  wings, 
lies  in  the  consideration  of  two  very  distinct  flow 
region*:  one  Is  a  nearly  planar  (2-0)  region  rait 
to  tne  wing  section  with  a  straamwise  length  scale 
comoarabla  to  the  typical  wing  chord  Co,  and  the 
otner  is  a  fully  3*0  domain  with  Its  size  compara¬ 
ble  to  the  wing  half-span, b.  They  are  the  "Inner" 
and  the  "outer"  region,  adopting  Van  Dyke's^) 
terminology.  Since  the  aspect  ratio  /ft,a*M^ls 
high,  the  flow  in  the  outer  region  sees  the  wing 
end  its  near  waka)  as  a  line  of  singularity  -  the 
lifting  line  -  in  the  leading  approximation. 

It  may  be  clearer  for  the  present  exposition 
to  begin  with  the  more  familiar  llftlng-llna  solu¬ 
tion  for  the  outer  region.  It  is  also  convenient 
to  write  dawn  here  the  definitions  of  the  dimension 
less  Inner  and  outer  variables.  For  the  outer 
region,  we  wi 1 1  use  the  variables 

*  f  V/Tmi  * ■  /*//*,  *«**. 

?  z  a  f/**  14 c,  . 


2.1  The  Induced  Upwash  Field 


The  perturbation  velocity  ootential  #, ( 7,  y, *  ) 
representing  the  combined  system  of  the  bound  and 
trailing  vortices  in  the  outer-flow  region,  satis¬ 
fying  the  Prandtl-Glauert  equation,  Is 
.» 


j  :  I  f  ^(7. ) —  r  t  *_i5l£L>1  if 

*.■  |C +(?-?.)***'  |* 


(2.4) 


The  higher  approximation  of  the  outer  solution 

U.5) 

must  account  for  not  only  the  finite  aspect-ratio 
affect  as  In  Van  Dyke's  "third-order"  theory^), 
but  also  the  nonlinearity}’**’®)  although  the  latter 
wj,ll  affect  little  the  3*0  corrections  to  the  sur¬ 
face  pressure  of  interest.  For  the  type  of  center 
lines  considered  below,  wa  have  either  IsAjf  or 

Xj  *171  .  matarnA/fmZ  • 

Wa  point  out  that  the  roll-up  tendency  of  the 
tral 1  log- vortex  (TV)  sheet  far  behind  the  wing  is 
controlled  by  the  trailing  vort I c I ty 4t,"V/%y.  Thus, 
for>ft,)Pf,  tha  assumption  of  a  flat  TV  sheet  im¬ 
plicit  in  Eq.  (2.3)  is  seen  to  be  asymptotically 
correct,  as  long  at  I  ■  0(1),  irrespective  of  the 
degree  of  nonlinearity  in  the  inner  flow  region 
(next  to  the  wing  section).  Therefore,  this  assump¬ 
tion  will  not  preclude  the  proper  analysis  of  the 
large  deflection  of  the  TV  sheet  in  the  x'-z'  plane 
next  to  the  wing  section  in  a  ful 1 -potent ia I  formu¬ 
lation.  (57) 


and  for  the  inner  region,  we  use  the  variables 

***'/*. 


0  s  ti/dPu.c,. 


(2.2) 


In  the  aoove,  (x.y.z)  is  the  right-handed  Cartesian 
coordinate  system  with  the  x-axis  pointing  in  the 
free  stream  direction  and  z-axls  in  the  lift  direc¬ 
tion;  whereas  (x’.y'.z1)  is  a  right-handed  ortnogo- 
net  curvilinear  coordinate  system  with  z'  •  z  and 
y'  *  0  being  the  center  line  of  the  wing  planform 
*  •  «c(y) .  and  x'  Is  the  distance  from  the  center 
line  measured  In  x-y  plane  fcf^_  Fig.  2).  Itjs 
assumed  chat  the  reference  center  line  Is  located  ” 
completely  in  the  x-y  plane,  l.e.  we  are  consider¬ 
ing  a  planar-wing  problem.  Note  that  Eds.  (2.1) 
are  in  effect  the  Frandtl -Glauert  variables,  with 
che  scale  for  the  perturoetion  potential  so  chosan 
:o  match  that  of  the  inner  solution;  and  that  Eqs. 
(2.2)  are  si.xply  the  set  suggested  by  the  transonic 
similarity  law  for  the  component  flow,  using  the 
naif  reference  chord  co/2  as  the  basic  langth  scale 
*nq  %  ■  H,co*A  as  the  component  free-t.tVeam  ve lor 
Itv.  were  0,  is  withal  being  taken  a 

typical  (absolute)  angle  of  attack,  or  the  wing 
sec:  on  chickness.  F or  all  practical  purposes,  wa 
assume  y  and  as  belong  to  the  same  order.  For 
convenience  in  the  subsequent  discussion  pn  the 
-elated  classical  work,  we  will  employ  a  dimension¬ 
less  inner  variable  i‘  as  an  alternative  to  the 
2  i n  Eq .  (2.2)  name  I y , 


A  ,  m*  • 

t  «  4*  /  c,  , 


The  most  crucial  knowledge  from  Eq.  (2.4)  is 
the  flow  behavior  in  the  vicinity  of  the  lifting 
line,  i.e.  2-Ee(9)^0  and  I-*0,  from  which  the 
Induced  velocity,  hence  the  incidence  correction, 
can  be  determined.  For  center  lines  Involving  only 
straight  segments,  the  behavior  in  question  can  be 
written  as 


Har$» 


am  .  .1  fA 


-7  £■ 

where  *  |Jf f)  .  and  £'*  a  ♦./'-W.1  i'  Note 

that  the  result  has  been  expressed  in  the  Inner 
variables  x,  y  and  z  to  facilitate  easy  identi¬ 
fication  for  subsequent  discussions  on  the  upwasn 
effect;  the  last  term  is  independent  of  x  and  z 
with 


*1  ^>>V^I  -Z.nA  ,  (2.6b) 

where  •  0  for  a  straight  oblique  wing  'Seft  y 
and  , 

A.mVcr*\-(  r  -l/l) 

*,(?  “  l  (f.  -f  )*'L 

-  Si*A  Sf»  (% .  f)  J  4f,  (2.6c) 


for  a  symmetric  swept  wing  3e*W).  in  above, 
*•"  A  r$MA/f  I  -  <*«  ,  ana  eel  A  »  cosA/Jl-at*  . 


12.3! 


Observe  Chet  f*Js  rfj \/4y  -  see  A  df-iUy  . 

The  first  terra  on  the  right  of  Eq.  (2.6a)  gives 
the  concentrated  vortex  representing  the  wing  sec¬ 
tion,  the  second  group  of  terras  proportional  to 
*7,/*?,  absent  in  the 

classical  analyses,  results  from  the  non-vanishing 
spanwise  component  of  the  near-wake  vorticity. 

This  component  of  the  vorticity  induces  a  logar I th- 
mically  infinite  upwash  at  the  center  line,  as  may 
be  anticipated  from  a  consideration  of  the  local 
irrotationa)  2-D  flow  (cf.  Fig.  3).  This  is  re¬ 
flected  in  the  logarithmic  dependence  of  2  on 

the  aspect  ratio  .  The  terms  in  Eq.  (2.6a) _ 

which  can  be  identified  with  the  proper  incidence 
correction  (to  the  inner  solution,  after  properly 
matching  the  inner  and  the  outer  solution)  are 


7) 


i.'  (z+ze .  -i  i5  i**)  <2- 

[The  terra  with  the  factor  (In  2)  will  not  be  pre¬ 
sent  if  the  basic  length  scales  for  x',  etc.  is 
taken  to  be  co  instead  of  £co).  The  logarithmic 
terra  in  gives  a  significant  contribution  of 

the  near  wake  which  Is,  in  turn  controlled  by  the 
yaw.  Similar  (log  4,)  terms  will  rise  from  the 
center  line  curvaturel1**)  (or  from  a  time-depend¬ 
ent  7,  _  _)  • 

Two  observations  are  essential  in  subsequent 
discussions.  First,  the  part  of  the  induced  upwash 
Held  shown  in  Eq.  (2.6a),  which  varies  with  £  and 
z,  consists  of  the  concentrated  vortex  as  well  as 
the  logarithmic  behavior  of  the  near  wake;  second, 
since  log£,  results  from  the  logarithmic  singu¬ 
lar!  ty  in  the  upwash,  the  log/R,  term  would  have 
been  absent  from  the  ir  Eq.  (2.6b)  and  expres¬ 
sion  (2.7)  above,  if  the  induced  incidence  correc¬ 
tion  were  taken  simply  (and  incorrectly)  as  the 
finite  part  of  the  upwash  field.  In  passing, 
sTmSTar  results  have  also  been  derived  for  a 
straight  obi ‘cue  wing  by  Cook'18), 
transformation  techniques. 

2.2  Related  Classical  Analyses 


using  he  1 1  in 


Among  the  earlier  methods  applying  the  lift¬ 
ing-line  idea  to  swept  wings  (see  Ref.  3,  28-30  for 
a  thorough  revtewl,  the  most  well  known  is  perhaps 
that  of  Wei  singer*-”).  As  pointed  out  by  Jones  and 
Cohen (3),  Wei  singer's  method  does  not  recover 
correctly  the  limit  for  the  infinite  aspect  ratio. 
We  may  recall  that  in  Weisinger's  method,  the 
bound  vortices  at  each  wing  section  are  represented 
by  a  concentrated  vortex  at  the  1/4  chord  point  and 
the  induced  downwesh  (negative  upwash)  Is  computed 
at  the  3/*  chord  point.  Nevertheless, the  logarith¬ 
mic  dependence  on  mentioned  should  be  found 
also  in  solutions  by  the  Weisinger  method  (applied 
to  high-aspect-ratio  wings),  since  an  upwash 
formula  equivalent  to  Eq.  (2.4)  was  used  therein. 

On  the  other  hand,  the  use  of  the  3/4  chord  as  an 
upwash  control  point  requires  a  uniform  upwash,  or 
one  varying  linearly  in  x' ;  but  this  requirement"' 
is  not  met  by  Eq.  (2.6a).  In  addition,  its  valid¬ 
ity  also  imply  the  assumption  of  a  local  2-0 
component  flow  devoid  of  the  near  wake,  which  Is 
again  contrary  to  Eq.  (2.6e).  I ts  fei lure  in  the 
high-aspect-ratio  range,  therefore,  should  not  be 
too  surprising. 

This  comment  does  not  apply  directly  to  the 
vortex- lattice  method^'  ,32)*mploy ing  tne  3/4  chord 
as  a  control  point  for  each  lattice  (panel).  This 
it  because  the  large  upwash  and  spanwise-vorticity 


jump  associated  with  the  logarithmic  singularity 
-  -aepend  mainly  on  the  vortex-shedding  rate  of  tne 
— hndfvidual  lattice,  whose  contribution  is  made 
smaller  as  the  lattice  number  increases.  [Error  cf 
such  a  nature  will  nevertheless  remain.) 

Solutions  to  elliptic  lifting  surface  at  yaw 
have  been  given  early  by  Krienes  in  Ref.  33,  wnere 
solutions  to  the  inverse  lifting  problem  were 
superimposed  to  determine  lift  and  moments  of  a  5:’ 
elliptic  flat  plate  at  yew  (a  helpful  delineation 
of  Krienes'  analysis  in  the  unyawed  case  can  be 
found  in  Ref.  3  i.  Of  the  five-term  truncated 
_  series  used  therein,  three  terms  were  syimnetric 
spanwise,  it  is  not  clear  if  the  remaining  terns 
could  adequately  describe  the  asymmetric  span  load 
of  interest.  Generalization  of  the  lifting-line 
theory  for  a  wing  with  a  curved  center  line,  as  well 
as  a  wing  1^  side  slip,  has  been  considered  by 
Dorodni tsyn '34)  who  noted  the  significance  of  the 
logarithmically  large  upwash  effect  due  to  yaw.'*' 
However,  the  analysis  involved  the.3/4  chord  con¬ 
trolled  point,  as  in  Weisinger's '2°) ,  and  therefore 
the  3*0  effect  analyzed  cannot  be  expected  tc  be 
asymptotical ly  correct.  The  results  were,  in  any 
case,  restricted  to  small  departure  from  a  straight, 
unswept  center  line  in  the  linear  regime. 

A  case  of  curved  center  line  has  been  studied 
by  Thurber'35)  who  considered  a  crescent-moon-snaped 
wing  from  the  viewpoint  of  an  asymptotic  theory; 
however,  the  solutions  for  the  region  around  the 
wing  section  and  the  matching  problem  were  not  con¬ 
sidered  (the  upwash  calculation  also  contained 
errors),.  Oscillating  high-aspect-ratio  wings  with 
curved  center  lines  have  been  treated  by  the  senio- 
author  in  Ref.  36,  where  the  basic  reasons  for 
the  appmarance  of  a  logarithmically  large  upwash  are 
explained. 

2.3  The  Asymptotic  Approach  and  the  Basic  ~eatu-es 
of  the  Theory  ~ 

The  works  under  review  adopt  the  basic  ap¬ 
proach  initiated  by  Van  Dyke (27)  who  considered 
Prandtl's  lifting-line  theory  as  the  leading 
approximation  in  an  asymptotic  analysis  for  a  hign- 
aspect-ratlo  wing.  While  this  more  formal  approach 
may  break  down  near  the  tip  of  a  rectangular  or 
taper  wing  (where  the  asymptotic  expansion  is  non¬ 
un  iform,as  is  quite  well  known),  the  methods  allow  a 
systematic  improvement  by  successive  approximations, 
and  a  clearer  delineation  of  the  various  competing 
3-0  effects.  Thus  Van  Oyke  is  able  to  extend 
Prandtl's  work  for  a  straight,  unswept  wing  to  a 
higher  order'’8'.  On  account  of  the  tip  nonunifor- 
ml ty ,  however .  the  third-order  results  of  Ref.  37 
does, not  predict  correctly  the  total  lift  and  crag; 
there  is,  in  addition,  an  error  related  to  the  Kutta 
condition,  to  be  discussed  later. 

An  example  of  the  asymptotic  analysis  of  this 
type,  representing  a  significant  departure  from  the 
assumption  of  a  straight,  unswept  center  line,  is 
given  in  Refs.  12  (a)  and  (b)  for  an  oblique  wing 
in  the  framework  of  the  linear  theory  of  an  incom¬ 
pressible  flow.  This  linear  analysis  snares  the 
same  outer  solution  and  the  corresponding  upwash 
distribution  in  the  subsequent  transonic  oblique¬ 
wing  *na lyses*' 4) _ _ The  rather  encouraging 


This  work  was  called  to~our  attention  bv  Ors.  G. 
Daforno  and  R.E.  Melnik  of  the  Grumman  Aerospace 
Corporation. 
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comparison  with  the  more  exact  wing-panel  method 
given  therein  suggests  that  equal  success  may  be 
obtained  for  the  transonic  swept-wing  analy-~~ 
ses('3-18)  at  least  for  the  oblique  wings*. 

(!)  Departure  from  the  classical  analysis 

Aside  from  the  interesting  dependence  of  the 
induced  upwash  on  the  span  loading  and  on 

the  sweep  angle  A  ,  to  be  delineated  more  fully 
later,  there  are  other  features  in  the  formulation, 
which  represent  important  departures  from  the 
classical  theory.  One  of  these  is  the  need  for 
considering  the  spanwise  component  of  the  trailing 
vorticity  between  the  trailing  edge  and  the  down¬ 
stream  infinity  in  the  reduced  inner  problem.  Such 
a  need  is  made  evident  by  the  flow  behavior  shown 
in_Eq.  (2.6a)  where  the  second  term  gives  a  Jump  In 
across  the  wing  trace,  being  proportional 
to  -  taoA/Rr'dPv/f/ •  This  means  that  the  local 
component  flow  cannot  be  determined  by  a  simple  in¬ 
cidence  correction  to  a  wakeless  2-D  component 
flow  as  in  Prandti's  original  analysis.  The  solu¬ 
tion  presented  for  the  straight  oblique  wing  in 
Refs.  12  (a)  and  (b)  confirms  the  importance  for 
treating  the  vorticl ty .iump.at  the  wing  trace.  In 
the  more  general  cases'1 4 the  partial  differ¬ 
ential  equation  (POE)  governing  the  inner  region 
must  also  be  corrected  for  the  center- line  curva¬ 
ture  and/or  for  a  3*0  compress ibi I i ty  effects  (see 
below) . 

One  must  observe  that,  as  the  component  Mach 
number  m,  approaches  one  (from  below)  ,  the  sweep 
angle  A  must  necessarily  decrease  toward  zero 
(unless  I ) .  According  to  Eq.  (2.6a).  a  smal  ler 
A  will  not  diminish  the  importance  of  che 
logarithmic  upwash  field  associated  with  the  span- 
wise  near-wake  vorticity,  but  reduce  at  the  same 
time  the  spanwise  component  of  the  wake-vorticity 
itself.  For  this  reason,  the  velocity  wi 1 1  remain 
essentially  continuous  across  the  wing  trace  in  the 
inner  problem  for  transonic  swept  wings. 

(i I )  Kutta  condition 

The  Kutta  condition  at  the  trailing  edge  plays 
an  imoortant  role  in  the  present  theory.  While  an 
infinite  (.integrable)  singularity  in  the  3-0  cor¬ 
rection  to  the  surface  speed  could  suggest  only  a 
local  breakdown  (nonuniform)  of  the  asymptotic 
solution  near  the  trailing  edge,  it  gives  an 
erroneous  correction  for  the  entire  wing  section. 
This  is  because  the  latter  represents  an  unwanted 
eigen  solution  to  the  homogeneous  boundary-value 
problem,  as  is  shown  quite  clearly  by  the  explicit 
inner  solution  in  Refs.  12  for  an  oblique-wing. 
Similar  problems  appear  for  the  unyawed  straight 
wing  at  the"tnird  order’where  the  particular  solu¬ 
tion  to  the  inner  problem  given  originally  in  Refs. 
27  and  37  violates  the  Kutta  condition, as  being 
uncovered  and  corrected  by  KIda  and  Mlyai'3«). 

2.4  Behavior  of  the  Induced  Upwash  on  Sweet  Wlnas 

_ The  pert  of  Induced  upwash  distribution  in  the 


The  scales  in  the  span  loadings  shown  in  Fig.  3 
of  Ref.  12(a)  and  in  Figs.  A-6  of  Ref.  12(b) 
are  incorrect;  they  should  all  oe  reduced  by 
a  ‘actor  of  10. 


form  of  2  is  dominated  by  the  term  ~(yt rr'(i-M„‘)‘.* 
•Sin A.  d  ff/ay  In  (J,  -  m.1  M,)  which  wi  1  (.general  ly 

reduce  the  downwesh  on  an  aft  (or  swept-back  _ _ 

wing)  panel  where  d/J/Vy  <  Q  and  augment  the 
downwesh  on' the  forward  panel  where  >  o .  For 

most  span  loadings  of  interest,  Jf,/dy  tends  to 
infinity  at  the  tip  J  •  1  I  like  (l-y*  ),  this, 
together  with  the  logarithm  1  nyo  I  v  InqJ^  S 
and  (i- }*iJeads  to  a  maximum  of  0(  /H"**)  near  the 

tip  of  a  swept-back  panel.  _ 

there  is  also  a  minimum  in  3  hear  the 
tip  on  an  swept-forward  panel  (where  the  dawnwash 
reaches  a  magnitude  of  OfA?**4)).  Although  the 
magnitude  of  the  Induced  velocity  becomes  infinite 
at  the  tip,  2  reverses  Its  sign  and  vanishes  at  a 
span  station  extremely  close  to  the  tips.  This 
tends  to  provide  a  reasonable  description  of  the 
span  loading.  In  spite  of  the  local  breakdown  at  the 
tip.  The  latter's  region  of  nonuniformity  is  esti¬ 
mated  to  be  y  t  i  s  OCX^UX,  ). 

The  asymmetrical  span  loading  contributed  by 
terms  proportional  to  -  J f$/Sy  in  £  are  respon- 
ible  for  the  unbalanced  rolling  moment  of  an 
oblique  (pivoted)  wing  ,  unless  twist, wing  bend, or 
special  pivot  location  is  introduced  in  the  de¬ 
sign^  II).  The  extent  to  which  this  asyntnetry 
depends  on  the  yaw  angle  A  and  on  the  type  of  the 
basic  span  loading  ftc/)  is  illustrated  in  Fig.  k 
where  the  upwash  functions  £ (y)  for  obi Igue  wings 
with  an  elliptic  and  with  an  extended-span  (ES) 
distribution  in  ^(9)  are  shown  for  *o,  22.5 * 
and  1*5*.  These  g’s  were  computed  from  Eq.  (2.6) 
with  2(7)  normalized  by  the  mid-span  value  of  fm, 
and  the  (reduced)  aspect  ratio  taken  to  be  such  that 
“  8.1*.  We  note  that  the  graphs  can  be 
corrected  for  other  aspect  ratio  by  simply  adding 
to 


/  5mA. 


rj(  f)-ln 


A,  stcAl 

B.4  I  . 
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The  ES  load  considered  has  a  root-bending  moment 
equal  to  that  of  an  elliptic  load  for  the  same 
lift  but  lesser  drag'1*0'.  The  span  however  Is 
longer,  hence  the  name  "extended  span".  For  the 
ES  distribution  shown,_ the  extended  span  is  1.15 
times  that  of  the  elliptic  one. 

We  note  in  passing  that  at  A.  “  0,  the  ES 
load  induces  a  lineal  distribution  in  the  upwash, 
equivalent  to  the  effect  of  a  wash- i n ,  and  the 
asynsnetry  caused  by  the  yaw  wilt  further  increase 
this  effect  on  the  downstream  side,  turning  a 
regular  downwash  into  a  (positive)  upwash  on  the 
aft-panel  tip.  The  amplification  of  the  induced 
velocity  in  the  out  board  region  of  each  panel,  and 
the  existence  of  a  maximum  and  a  minimum  in  the 
manner  discussed  earlier,  are  quite  apparent  from 
the  figure. 

The  induced  flow-angle  correction  will  depend 
also  on  the  center- line  geometry.  Fig.  5  shows  the 
span  distribution  of  normalized  2  *  £e0  for  sym¬ 
metric  swept-foward  and  swept-back  wings,  at 
different  degrees  of  sweep.  It  is  assumed  in  the 
calculations  that  the  span  loading  P0{y)  (based  on 
the  strip  theory)  lselliptic,  and  the  (reduced) 
aspect  ratio  4?,  s  / 1 - •»*  at,  is  9.26.  Results  for 
five  values  of  A  (0,  *  22.5*,i<*5*)  are  shown.  At 
M  —  ■  0.755,  A  *  t  22.5*  and  1 1*5*  corresponds  to 
•  tl5.2*  and  *  33-33*.  respectively. 

The  ‘ami  liar  fact  aoout  the  need  ‘or  a  wash¬ 
out  on  a  sweot-becx  wing  is  quite  evident  from  the 
large  reduction  in  |  £  *  2 J |  near  the  tip  shown. 


it 


The  usefulness  of  adopting  a  forward  sweeo  ( A  <  0} 
in  counteracting  the  above  effect  is  also  quite 
evident  from  Fig.  5. 

For  value  of  4 C,  other  than  9-26,  the  cor¬ 

rection  to  2  indicated  in  Eq.  (2.19)  again 
applies.  The  behavior  near  the  tip  of  a  swept  wing 
is  again  similar  to  that  observed  earlier  for  the 
oblique  wing;  but  in  approach  the  center  corre¬ 
sponding  to  the  wing  appex,  the  magnitude  of  the 
induced  upwash  becomes  unbounded.  This  divergence 
is  expected  from  2*  which,  upon  9  -*  0,  yields 


This  is  well  borne  out  in  the  k5®  swept-back  case, 
for  which  the  downwash  near  the  center  is  more 
severe  than  other  cases  shown.  Unless  JJ(0)  •  o, 
the  theory  will  break  down  at  y  ■  0.  Comparing 
{Z  +  R'.)fiC  with  the  leading  term  in  Eq.  (2.6a), 
in  conjunction  with  Eq.  (2.9),  reveals  readily 
that  a  complete  break  down  occurs  where  (t-  Mi  )'Vl  • 

•  A?' jin  A  f$(? )  /j  becomes  of  unit  order,  i.e. 
at  '  _ _ 

r  SinA  ?[(o)  1 

7  =  0 - — 1 -  (2.10) 

Lo-O*. 

This  nonuniformity  will  be  discussed  further  later 
in  the  context  of  the  nonlinear  transonic  flow. 

In  passing,  we  note  that  unlike  the  upwash  on  a 
swept-back  wing,. the  S  •  Si  near  the  *P*X 
station  of  a  swept-forward  wing  does  not  become 
excessively  large  until  y  comes  very  close  to  the 
origin.  Thus,  breax  down  of  the  theory  near  y  »  0 
may  be  less  severe  for  the  swept-forward  case. 

3.  TRANSONIC  SWEPT  WING:  THE  INNER  PROBLEM 

We  consider  high  subsonic  flight  near  the _ 

speed  of  sound ^  As  is  well  known,  the  performance 
ofli  transonic  wing  is  rimfted  by  the  rapid  "drag 
rise  as  H,^»l .  With  a  moderate  sweep,  however  uone 
can~6r frig  the  f Tight  Mach  number  of  a  high' aspect 
ratio  wing  closer  to  one,  while  keeping  the  compo¬ 
nent  flow  around  the  wing  section  at  the  verge  of 
the  drag  rise. 


3. 1  Basic  Parameters  and  POE 

The  basic  component  flow  in  this  situation  1$ 
necessarily  non  I  inear  and  of  a  mixed  (hyperbolic- 
elliptic)  type;  if  the  wing  section  Is  sufficiently 
thin,  this  component  flow  at  each  span  station  wi.lj 
obviously  be  control  led 'by  "the  transonic  similarity 
parameter  based  on  the  component  Mach  number 
M„  -  f^dotA 


is  that  which  will  keep  K„  at  the  unit  order,  while 
does  not  exceed  one.  Let 

©*4.  ^.|b> 

A 

The  foregoing  condition  is  then  equivalent  to  the 
requirement  ©»0(1).  Finally,  there  is  a  reduced 
aspect  ratio  controlling  the  3*0  influence,  the 
reciprocal  of  which  is 


Note  that  et  M,  *  (H,,K,,)  ^  .  The  asymptotic 

analyses  of  Refs.  13*18,  including  the  unswept  case, 
pertain  to  the  limit  €-w0,  while  keeping  K„  and 
0  (also  a™  e"1  )  fixed.  With  the  additional  re¬ 

quirement  M— <  1,  the  conditions  on0  ana  K„ 
amount  to 

©*  <  K,  *  0(1 )  (3.2) 

When  the  center-line  curvature  is  taken  to  be 
zero  locally,  and  the  disturbance  is  assumed  to  be 
sufficiently  weak,  the  PDE  in  the  flow  region  next 
to  the  wing  under  conditions  just  mentioned  can  be 
obtained  from  an  expanded  form  of  the  full  potential 
equation  in  the  rotated  coordinates  (x',y’,z:)  as 

0-M?)*,y  *  w»A  ky  - 

-  *i-  «■•••  =  0,  (5-3) 

where  the  terms  omitted  (•••)  can  be  shown  to  be 
relatively  small  on  the  basis  of  the  transonic 
smal  1  -disturbance  assumption  (i.e.  I-M,*  *  O'*"*) , 
p  *  OtU.  *.*"*)  ,  »/**■  ■  *»*'  *  i  •  o (*■'•)]. 

Except  for  the  third  term,  the  equations  shown  above 
are  familiar  from  the  transonic  smal 1 -di stur- 
bance  (TSD)  theory (J,0-«4) .  The  added  term, 

a»A  results  from  a  3*D  compress i bi  I  i ty 

correction  to  the  velocity  divergence  (owing  to  the 
spanwise  variation  of  the  density)  and  will  re¬ 
main  also  in  the  linear  regime,  unless  M„  or 
tanyi  vanishes.  In  terms  of  the  inner  variables 
[cf.  Eq.  (2.2)]  the  above  equation  can  be  written  as 

h  (Mf-  ¥  #]  *  +  -  '3  •  U) 

with  the  remainder  being  comparable  to  €*  .  Thus 
the  3*0  compressibility  correction  appears  as  a 
correction  to  the  Von  Karmdn  equation^0)  and  the 
reason  for  introducing  the  reduced  sweep  angle 
0a  A/  ecv*  is  made  apparent.  A  pressure 
coefficient  may  be  computed  from  the  solution  as 


familiar  from  the  transonic  smal 1 -disturbance 
theory  (l,0*‘,3) ,  on  account  of  the  stronger  spatial 
attenuation  for  disturbances  in  a  3*0  flow  field, 
it  is  possible  to  match  this  predominately  2-D 
non  I  inear  •  inner  flow  with  a  basically  J_lnea£  outer 

flow  described  earlier.  As  observed  i^TT  the _ 

swep.t_a.ngle  A  must  be  confined  to  a  relatively  nar- 
ow  range,  because,  if  A  is  too  small,  M„  will 
exceed  the  drag-rise  value  of  the  wing  section,  and 
if  A  I*  too  large,  the  M  would  become  so  far  from 
unity  that  the  advantage  of  a  high  pressure  co¬ 
efficient  will  be  lost.  The  A  range  of  interest 


t,  vy  -n 
ices  A  a,  <pt  i- 


With  a  remainder  comparable  to  e* .  Note  that 

as  long  as  0  %0(l) ,  i.e.  A-0(«1'1),  Cp  is 
unaffected  by  p?  . 

We  note  that  in  the  transonic  similarity 
parameter  defined  in  Eq.  (3-la),  the  denominator  is 
taken  as  the  first  power  of  M„  as  in  Murman  and 
Cole '*3)  also  Murman  and  Krupp'  .  This  choice 
gives  a  definite  improvement  in  the  accuracy  of  the 
critical  flow  speed,  hence, the  wave  pattern  near 
the  sonic  boundary,  as  pointed  out  in  Ref.  Ik  (b) 
and  in  Ref.  k6.  One  also  may  note  that  most  of  the 


"  cot  A.  "  factors  appearing  in  Eqs.  (2.2),  (3.4)  * 
(3-7),  except  that  in  (1-M,?  )  of  K,  ,  can  be 
replaced  by  unity  because  of  ©  »  0(1).  However, 
its  retention  in  Eqs.  (2.2)  and  (3-5)  will  allow  a 
greater  range  of  the  sweep  angle  \  for  the  valid¬ 
ity  of  the  theory  l'4b) . 

3.2  Conditions  on  Wing  Surface.  Wino  Trace  and 
ShocixT 

for  an  impermeable  thin  wing  surface  with  the 
leading  and  trailing  edges  in  the  wing  plane  de¬ 
scribed  by  x  ■  S(y)  and  x  «  6(y),  respectively,  the 
ordinate  can  be  written  as 

*'  =  4[<*£V*,y)  .  . 


If*)] 


A  ♦  ^  A 

where  Z  ,  Z$  ,  and  I  are  all  of  unit  order,  and 
the  in  signifies  the  two  different 

distributions  for  the  upper  gnd  the  lower  surfaces. 
The  functions  2,  (y)  and  I  (y)  represent, 
respectively,  a  wing  twist  and  an  (upward)  wing- 
bend;  they  are  introduced  to  allow  control  of  the 
3-0  effect.  The  wing  boundary  condition  can  be 
transferred  to  the  wing  plane  (z  «  0)  to  read 


where  the  subscript,"*"  refers  to  the  wing.  On 
the  wing  trace  at  *  <•  0  corresponding  to  the 

vortex  sheet,  the  sma I (-disturbance  approximation 
gives  the  requirements  that  pressure  jump  and  the 
normal  velocity  are  zero.  In  view  of  Eq.  (3-5), 
this  leads  to  the  condition  that  the  two  components 
of  the  velocity  perturbation,  $1  and  are 

required  to  be  continuous  across  the  wing  trace. 
Therefore,  in  solving  for  the  3*0  correction  for 
the  inner  solution, we  can  treat  the  problem  as  be- 
ing  wakeless.  This  simplicity  follows  from  the 
fact  already  noted  that  the  spanwise  component  of 
the  wake-vort ic i ty  is  not  large  enough  to  warrant 
consideration  in  as  much  as  is  small. 

The  PDE  (3 - **)  is  of  a  mixed  type  depending  on 
K „  and  iff  ;  the  characteristic  surface  of  the 
hyperbolic  region,  x  «  xc(£,y)  is  given  by(‘*7) 

K,-**,;  §f=  -(§1-7.  <3.8) 

in  the  limit  •-♦0,  we  identify  the  critical  speed 
sor  the  2-0  component  flow 

*  *  yr/. 

The  POE  (3-4)  admits  weak  solutions  with  disconti¬ 
nuity  surface  consistent  with  the  Rankine-Hugon iot 
shock  relations  under  the  TS0  approximation.  One 
of  toe  jump  condition  is  of  the  same  form  as  the 
cnaracter i st ic  equation  (3.3),  with  x*t/l  there  re¬ 
placed  by  x  (the  superscript  "D"  refers  to 

the  shock  discontinui tv)  and£j-  replaced  by  the 
ar'  thmecical  -seen  of  across  the  shock.  The 
otner  cond i : ion  governing  the  shock, corresponding 
to  the  continuity  of  the  tangential  velocity, can 
oe  simply  taken  to  be  the  continuity  of  I. 


As  seen  above,  the  3*0  correction  considered 
represents  a  small  peturbation  of  the  basic  2-0 
component  flow.  The  analysis  can  be  mede  more  sys¬ 
tematic  by  solving  for  tne  coefficients  in  the  _ 
asymptotic  expansion  for  small  g  : 


A  £ 

< p  -  Pt  '  *  9, 


(3.10a) 


with  a  corresponding  expansion  for  the  unknown 
shock  boundary 


(3.10b) 


where  a  weak  logarithmic  dependence  ofplmP  X,on«  is 
expected  (in  view  of  the  nature  of  the  upwesh  ex¬ 
amined  in  $2).  The  resulting  equation  of  particular 
interest  Is  POE  governing 

{[kw-m-)  &?]  &*&,- «">  iffh}  $, 

-  (3  - 1 1 ) 

The  nonhomogeneous  term  on  the  right  results  from 
the  3-0  compressibility  correction,  it  confirms  the 
remark  given  eerl ier  2) that  the  upwash  correc¬ 
tion  alone  does  not  give  the  complete  answer  to  the 
3*0  potential  problem. 

As  in  most  standard  perturbation  analysis  invol¬ 
ving  shock  waves,  the  jump  conditions  are  ana¬ 
lytically  transferred  from  x  »  x  D(z;y)  to  the  un¬ 
perturbed  boundary  x  «  x  § (2 •, y)  where  the  resulting 
condition  for  I  and  x?  read 


2®  -  <**•)  < 


■»*  ■»£  > 


(3. 12a) 
(3.12b) 


where^  ^  and[£  JJstands  for  the  arithmetical  mean 
and  the  difference  across  the  discontinuity  surface, 
respectively. 

3.4  The  Far  Field  of  the  inner  Solution 

At  a  point  sufficiently  far  from  the  wing  sec¬ 
tion,  the  solutions  and  ft  admits  a  behavior, 
apparent  from  the  governing  POE, 

%  ~  Tn  [>"  {ff  )  +  T  %**]♦-  <3  • ■ 1 

A  A 

^  2  d  e  r  a  |  ^  L  CT  A  a  f  *4 

Wher*  <3- '3b; 

As  well  known,  the  behavior  shown  in  (3.13a)  repre¬ 
sents  a  concentrated  vortex  with  a  circulation 
equal  to  the  potential  jump  at  the  trailing  edge. 
The  leading  term  in  Eq.  (3 . 1 3b )  's  simply  the  con¬ 
tribution  of -x©  Jkf-on  the  RHS  Eq.  (3-113  ,  using 
Eq.  (3.13a).  The  next  two  terms  in  Eq.  (3.13b)  are 
simply  the  two  components  of  the  induced  velocity 
to  be  determined  by  matching  with  the  outer  solu¬ 
tion. 

The  remainder  of  (...)  in  (3.1 3* )  Is  comoar- 
able  in  magnitude  to  the  doublet,  i  .e.  to  l£j''  , 
including  term*  proportional  to  and 

44* a*/  1 A | ♦  ,  as  is  quite  well  known  from  the  2-3 

TSD  analyses.  The  remainder  (...)  co  r  p ,  in  ( 3  •  3  s) 
is  comparable  to  unitv,  (to  be  more  precise,  to 
and  unity),  including  terms 
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proportional  to  rt  as  well  as  the  douDlet  strength 
in  the  far-fielc  o*  0O.  It  must  be  pointed  out 
that  these  remainders  are  not  essentia J_  for  formu¬ 
lating  the  reduced  problems  for  and  0|  ,  but  its 
use  at  Che  far  boundary  in  a  numerical  computation 
proves  to  be  very  helpful  ('3*  '*•'  . 

3-5  A  Line-Source  Effect 

Amonq  the  ioqari  thmic  terms  mentioned. is  a 
(pure)  line-source  term  its  existanee 

could  have  been  anticipated  from  the  appearance  of 
the  source-like  3*0  compressibility  correction  on 
the  RHS  of  POE  (3.4)  or  (3.11).  Ref.  14(a)  shows 
that,  for  a  straight  center  line,  the  source  strength 
5,(9)  can  be  explicitly  computed  from  one  of  the 
doublet  strength  belonging  to  the  far-field  expres¬ 
sion  of  S0;  the  result  is  also  applicable  in  cases 
involving  piecewise  straight  segments. t  This  and 
another  term  of  the  order  Jog  1^1  are  not  found 
however  in  Cook's  analys is ' '8) • 

The  existence  of  a  nonvanishing  Q,  ca^n  be 
readily  demonstrated  for  cases  in  which  can  De 

neglected.  Note  Eq.  .(3.11)  becomes 


( y  ,  i1  \  * 

v is**  ' 


(3.14) 

K,  3t  3?  0 

in  tne  far-fieid  (|5|-»ao).  In  this  case,  the  area 
integra.l  of  the  RHS  of  Eq.  (3.14)  can  be  equated  to 
the  total  volume  flux;  therefore, 

%  K*  Jj  *2*1  °  7 

1®  \im  ^  2®r,  tis  (3.  15) 

where  is  the  far-field  doublet  strength  of  0O 
associated  with  the  airfoil  thickness.  It  is  thus 
seen  that  the  swept  or  oblique  wing  will  presert  a 
line-source  effect  in  the  outer  flow,  which  must  be 
taken  into  account  in  a  third-order  theory  (in  Van 
Dyke's  sensed,  37))t  even  in  the  linear  case. 

This  line-source  will  not  affect  the  lifting  part 
of  the  linear  problem,  of  course. 

3.6  Hatching  with  Outer  Solution 

In  the  inner  variables  x,y  and  z,  the  behavior 
of  the  outer  solution  near  the  centerline,  Eq. 
(2.6a),  can  be  expressed  asu<*bJ 


~  C££.,)  [  ttn'ivfc,  2)  *  5  *  3  * 


xtr 

■f 


+  «(f  ►z'.jtf  ” t3. 16) 

where  we  have  omitted  ali  factors  of  Mn,  such  as 
Mn*^*  In  tan~'(Mn'il  *  /SiT,  z  \  since  its  error  can 
affect  the  determination  of  J*0((y))  at  most  to  the 
order  ot  or  which  in  turn  can  influence 
the  inner  solution  through  the  upwash  correction  to 
the  order  e*  ,  at  the  most.  The  third  term  on  the 
RHS  of  Eq.  (3.16)  confirm  the  anticipation  that  the 
spanwise  vorticitv  jump  can  be  neglected,  which  is 
seen  to  be  proportional  to  m  2/3 »  ,  hence  compar¬ 


able  to 


Comparing  Eq.  (3-16)  with  the  corresponding 
far-field  behavior  of  the  inner  solution  (0_  +tp,+ 
,  and  noting  jhe  slight  difference  in  the _ 
definitions  of  l  and  0,  we- can  now.,  ident  i  fy„ 
f*0(y)  s  ^t(y))  with  cos!/3jl  JTo(y)  a:  IT,  (7), 
and  the  induced  velocity  correction  with  (J+JJ): 

*  £  al  so  dependson 


A, 

c  C  ’ 

A  ‘  ’  _  _  )  ?3. 17, 

r-  'v  1  -  t -  rc  -  U2.©dr„j 

With  these  identifications,  the  inner  and  outer  so¬ 
lution  shown  in  Eqs.  (3-13)  and  (3-16)  are  seen  to 
match  [for  all  terms  snown,  except  the  third  (higr.- 
er-order)  term  in  Eq.  (3-16)].  Ref.  14  also  con¬ 
firms  the  matching  to  a  higher  order  in  l£|"'  and  e 
between  (S0  *£?,),  and  (0o  *  £  0 , )  ,  in  which  the 
strengths  of  the  line  doublet,  the  vortex,  and  the 
line-source  for  the  outer  0,  solution  are  determined. 

/N  A  . 

a  With  <£ ,  and  now  determined  in  terms  of 
!rj,(y).„it  is  assumed  that  the  inner  solutions  for 
20  and  S,  ,  wnich  fulfill  the  Kutta  condition  and 
are  continuous  everywhere  except  across  the  wing 
and  in  the  wing  trace  (x  > a ,  :  *  i  0),  are  uniquely 
aetermined  by  their  PDE's,  together  with  their  re¬ 
spective  conditions  on  the  wing  and  on  the  wing  cr ax, 
and  the  far-fieid  conditions  Eq -(3-13)  • 

The  problem  of  swept  wing  is  thus  reduced  to 
solving  first  a  strictly  2-D^systera  involving  the 
nonlinear  TSD  Eq.  (3.4)  for  0O  with  its  riant-hand 
member  omitted,  and  U  therein  replaced  by  So,  and 
next,  solving  another  2-D  (linear)  PDEAsystem  with 
the  PDE  (3,11)  for  the  3~D  correction  0,  .  The 
spanwise  variable  y  appears  in  the  two  (uncoupled' 
PDE  system  only  as  a  parameter  through  the  wing  and 
the  far  boundary  conditions.  The  solutions  at  each 
span  station  can  be  determined  independently  of  an¬ 
other,  except  for  the  spanwise  distribution  of 
?fo ( y )  and  4  j  ^  ,  which  can  be  obtained  after 
has  been  determined  from  an  adequate  numoer  of 
span  stations 

n. 

3- 7  Superposition  of  Similarity  Solutions  for  B. 


^  Owing- to  the  linearity,  the  3-D  correction  to 
0,  can  be  decomposed  into  separate  parts.  There  is 
an  important  class  of  wing  surface  geometrv,  for 
which  each  of  these  separate  parts  (after  being 
suitably  scaled),  has  similarity  solutions  indepen¬ 
dent  of  y,  as  does  the  basic  solution  0O.  Thus, 
the  reduced  2-D  equation  system  in  this  case  can  ce 
solved  once  for  al  1  span  stations.  [This  interest¬ 
ing  solution  stucture  does  not  appear  to  be  recog¬ 
nized  bv  Cook  in  her  anal  vs  I s( 18) .] 

This  wing  class  requires  that  the  wing  section 
at  each  span  station  be  generated  from  a  single 
airfoil  profile,  but  the  local  chord  c(y)  may  vary; 
local  twist  and/or  wing  bend  for  compensating  tne 
3-D  effect  are  also  allowed,  as  in  Eq.  (3-6).  This 
geometry  can  be  more  specifically  written  as 

*'r£[*cc;; ZTy) 

where  c($)  3  c(y)/c0.  Implicit  in  Eq.  /3-1 
that  the  center-line  x  »  z  *  0  is  the  common 
straight  axis  for  the  similar  wing  sections  at  dif¬ 
ferent  span  stations;  the  percentage-chord  location 
of  this  straight  axis  can  be  arbitrarily  chosen, 
however . 


To  describe  the  flow  structure  in  auestion,  we 
shall  rescale  x  and  z  by  the  local  chord  c;v), 
and  introduce  the  variables 


'=  */£.  z  =  i/t  ,  y*ecst.  ?  .  ?=  */;  IMS' 


The  similarity  flow  structure  admissible  under  Eq. 


(3-18),  with  due  allowance  for  the  upwash  correc¬ 
tions,  is  represented  by 

p  =  $,(*,?)  *■  6Qc'?:(x'.i)  -  t 


i’f(r)  *  1C?)  -©V/>]  %(?.z) 

A'  ^  (3.20a) 

where  c'=  dc/jp  ,  and  0O,  9,  and  02  are  independent 
of  y.  The  corresponding  shock  geometry,  when  an 
imbedded  supercritical  flow  region  appears,  takes 
on  a  similar  form 

*  =  £(*)  *  €<Qc'x?(7)  v 


♦  €. 


[k.  Cf<r)  - 

—  ,  (3.20b) 

The  third  term  of  Eq.  (3 • 20a)  is  introduced  so  that 
the  uowash  effect  is  transferred  through  the  solu¬ 
tion  Co  the  wing  surface  as  part  of  the  inci- 
dence  correction  to  fr0,  whereas  the  second  term 
60  C /<>>/  accounts  for  the  3-0  compressibility 
correction  to  the  POE  without  altering  the  wing 
boundary  condition. 

A 

The  POE  for  9n  is  the  same  as  that  for  0g  with 
£  ai.d  z  replacing  x  and  z;  the  linear  POE  governing 
3i  obtained  from  Eq.  (3- II)  is 

{[K. -«->&]  £.  *  *~t 

(3-21) 


i”1*; 

*  <m  zygl1]  V\?ll, 


(3.23c) 


where  9  m  tan'V  »/.'R,2r)r  f  These  results  are 

completely^recovejable  from  the  corresponding  re¬ 
sults  for  0O  and  0,  In 'the  more  general  case,  after 
observing  jff  #  ft  *  ,  £'  =  6Z?- 

-(<*)<)(.  Km)  (ir./irT  In!  ,  <3,  =  SSfff,  .  *tc. 

Mote  in  particular  that,  unlike  0,  in  Eq.  (3.13b), 
the  upwash  contribution  £i  T  appears  in 
neither  7,  nor  fy.  and  that  if*  ,  ft ,  ft  ,  ,  0%  , 

ate.  are  pure  constants. 

The  upwash  effect  now  appears  as  an  incidence 
correction  $)  in  the  reduced  problem  and  is 

determined  from  . . . . 


JZi;  =  fR  cf  - 

^  ’  4  *  >  2yR.tr 


(3.24) 


where  the  right-hand  member  results  because 

loST^ox  *nd  =  '  The 

02  satisfies  the  homogeneous  part  of  Eq.  (3.21). 

On  the  wing  portion  of  the  x-axis  (**io  )  , 

'z’K/li  -  i  zVkx  ,  tf’/ii  s  O  ,4nd  *9L/»t  *  1  . 

The  jump  conditions  for  a  shock,  if  it  exists, 
must  be  derived  from  Eqs.  (3.12),  not  from  the_ 
jump  conditions  for  the  weak  solutions  to  Eq. 

(3.21)-  Since  if,,  and  $,«0f $' &,  £?  e  z* 

+  €ilKm£f*  I  Z'»l  2  ft  ,  the  conditions  in 
question  can  be  deduced  from  Eq.  (3-12),  to  be^_ 
applied  at  the  unperturbed  shock  boundary  x  »  x^(z) 
We  give  below  the  conditions  for  f, 


with  Jic,  C,k  evaluated  from  Eq.  (3.17)=.  The 
source  strength  can  be  determined  with  the. help  of 
Green's  Theorem;  for  the  system  with  similarity 
structure  considered,  Q.  can  be  explicitly  evaluaed 
fromCtt) 

where  S’  and  B  are  the  leading  and  trailing  edge  lo¬ 
cations  in  5T,  respectively. 

3.8  An  Unsteady  Analogy 

For  wings  with  sections  profile  not  being  gen¬ 
erated  from  a  single  shape,  the  greatly  simplified 
solution  procedure  based  on  the  similarity  struc¬ 
ture  discussed  in  §  3.7  is  not  applicable.  An  un¬ 
steady  analogy  exists,  however,  which  provides  an 
alternative, and  perhaps  numerically  more  effective, 
solution  procedure  for  the  inner  problem  in  the 
more  general  case. 


Let 


t  s  //«€>  ,  f*  pi-<{KmC,i,  *3.25a,b) 


the  POE  (3.4)  then  appears  as 


(3-22a> 

«0.  (3.22b)  *(****-¥*«)*  ***  s  2<p*e-  f3’26) 


The  far-field  behavior  of  ?0,  7,  and  can  be  de¬ 
veloped  'rom  their  governing  POEiand  will  be  given 
aelow  more  fully  than  Eq.  (3-13),  inasmuch  as  they 
are  explicitly  used  on  the  far  boundary  in  the 
confutation  work  discussed  in  §5.  Here,  for  large 

p  u7  *i.  /R,  z  , 


<b 


£^r]  ^V3-23*) 

|  &  Imifi  - 

.7.  ITT  b  1  S' 

W2  ■  r r  . . 

•f!1  T  1  1^1 »  J 

*  *4  *  i3-:3U 


This  is  a  familiar  form  in  the  analysis  of  unsteady 
transonic  sma 1 1 -disturbance  flow  in  the  nonlinear 
,-egime(48,  49),  The  wing  boundary  condition  Eq. 
(3.7)  now  reads 


-n  *  a 
!£\  *2* 
dS  dt 


r  6  (  r  *!Z,  C 


(3.27! 


to  be  applied  at  z  *i0,  alV)<  *  <  S ti)). 

Thus,  an  analogy  of  the  inner  problem  eor  3  with  an 
unsteady  2-0  transonic  airfoil  problem  exists,  in 
which,  since  6  0  is  small,  the  unsteady  motion  is 
confined  to  the  neighborhood  of  the  cuas i -steady 
limit.  The  crucial  input  to  such  an  unsteady  prob¬ 
lem  is  of  course,  the  incidence  correct  ion  €jl(m' 

•  if.tTfl)  which  is  a  functional  of  7^20;  tne  latter 
can  be  replaced  by  an  approximate  distribution 


3 


a 


from  a  solution  of  Eg.  (3.26),with  JIT.  (£  t  set 
egual  to  zero  in  Eq.  (3.27) -  With  the  ADI  algorithm 
(49),  this  approach  promises  an  effective  compu- 
tational  procedure,  especially  in  capturing  shocks 
and  in  avoiding  the  ambiguity  associated  with  the 
reexpansion  singularity  at  the  shock  root  as  well 
as  the  problem  with  a  very  weak  shock  (discussed  in 
$  *••3). 


In  passing,  we  may  rarer*  that,  if  a  2-D 
shock-free  supercritical  profile  were  to  be  used  in 
the  3*0  wing  design,  the  unswept  (straignt)  wing 
appears  to  posess  one  advantage,  (although  there 
are  other  obvious  drawbacks)  since  the  nonhomoge- 
neous  terms  of  the  PDE  (3.4)  vanishes  in  this  case, 
and  the  remaining  3*0  effect  can  be  eliminated  bv  a 
proper  twist  distribution 


Numerical  results  have  been  obtained  for 
oblique  wings  via  this  approach  in  an  unpublished 
work  by  T.  Evans  and  the  senior  author,  and.  compare., 
reasonably  well  with  those  obtained  from  the  local- 
similarity  solutions.  Some  of  Evans'  nonlifting 
examples  nave  been  discussed  in  Ref.  14  (b). 

4.  LIMITATIONS  AND  COMMENTS  ON  THE  THEORY 

The  asymptotic  nature  of  the  analysis  not  only 
restricts  the  application  mainly  to  high-aspect- 
ratio  wings,  but  also  lead  to  a  number  of  nonunifor 
mities,  i.e.  regions  of  local  breakdowns,  which 
limit  considerably  the  usefulness  of  the  theory. 
There  has  also  been  questions  raised  on  the  unique¬ 
ness  and  existence  of  solutions  to  the  reduced 
problems.  Comments  on  these  and  other  questions 
are  given  below. 

4. I  The  Uniqueness  Problem  and  Shock-Free  Super¬ 
critical  Airfoil  ...  .  _  _ 

The  existence  and  uniqueness  of  the  reduced 
problem  for  0O  is,  of. course,  tne  same  in  the  2-D 
TSO  theory '*•<),  50,  51 ) .  For  the  uniqueness  of  the 
3-D  correction  c?  ,  the  problem  may  be  considered 
being  equivalent  to  asking  whether  the  homogeneous 
part  of  PDE  (3.4)  has  a  solution  with 

”  0.  It  is  essential  to  point  out  that 

the  swept  and  unswept  wings,  as  we lias  the  2-D _  . 

airfoil,  share  the  same  unfq'ueness 
problem  mentioned,  since  the  3~D  influence  enters 
only  as  a  nonhomogeneous  term  in  their  equations. 
Thus  the  uniqueness  proof  given  by  Cook'*®!  for  the 
lifting-line  theory  of  the  straight,  unyawed  wings 
in  transonic  flow  should  also  be  applicable  to  the 
swept-wing  case.  [The  idea  in  Cook's  proof  runs 
parallel  to  Morawetz's  earlier  work  on  2-D  mixed- 
flow  In  the  hodograph  plane(52),  but  differs  from 
Morawetz's  in  the  use  of  the  physical  plane  and  in 
allowing  circulation.] 

Two  important  questions  arise  concerning  the 
analysis  in  which  0-  represents  a"shock-f ree  super¬ 
critical  (component)  flow:  (i)  can  a  shock-  free 
perturbation  et,  exist?  (ii)  How  should  the  shock 
be  treated  in  an^asymptot i c  theory  for  small  C  in 
this  case,  if  e  0,  is  not  shock-free?  Definitive 
answers  to  these  questions  are  not  ava  iffbJb  to 
date.  In  fact,  the  proof  given  by  Cook('®)  aims 
chiefly  at  proving  uniqueness  for  the  perturbation 
of  a  shock-free  supercritical  solution  3^,  with  the 
implicit  assumption  that  the  perturbation  solution 
(9,  will  exist? 


*  One  must  recall  that  the  significance  of 
Morewetz'.s  original  uniqueness  study  was  concerned 
with  the  question  of  existence  (53),  where  she 
found  that,  to  maintain  an  imbedded  shock- free 
supercritical  flow,  the  airfoil  prof i le  variation, 
say  *  ■  t -ft?) ,  cannot  be  arbitrarily  prescribed  on 


[Of  course,  this  condition  can  be  useful  also  to 
the  oblique  and  swept  wings,  since  tne  contribu¬ 
tion  from  the  other  correction,  i.e.  the  3-D  com¬ 
pressibility  effect,  is  relatively  small  ir,  magni¬ 
tude  in  many  cases.] 

4.2  Breakdown  Near  the  Leading  Edge 

The  breakdown  of  the  asymptotic  analysis  of 
$  3  can  be  separated  into  two  categories.  One 
arises  mainly  from  a  smal 1 -disturbance  approxima¬ 
tion  which  must  fail  in  the  vicinity  of  a  stagna¬ 
tion  point,  or  where  the  x-gradient  is  singulars 
the  other  category  associates  more  directly  with 
the  high-aspect-ratio  expansion,  in  which  tne  soan- 
wise  gradient  is  assumed  to  be  small  and  must  fail 
where  Ts  singular.  The  nonuni  formi  tv  cf 

the  second  category  will  be  examined  in  $4.3. 

Stagnation  points  occur  near  both  leading  and 
trailing  edges,  the  problem  is  less  serious  witn  a 
sharp  trailing  edge  since  with  the  kutta  condition, 
the  region  of  this  breakdown  is  exponentially  small, 
as  Is  well  known,  for  an  airfoil  with  a  leading- 
edge  radius  comparable  to  t’C,  ,  2*is  proportional 

te  near  the  lcadingjedge  and  the  TSO 

theory'53,  54)  gjves  an  inf  ini  te  like 

(x-a)”,'?.  The  same  behavior  remains,  in  fact, 
even  if  the  cambe r,o r  absolute  incidence, is  comp¬ 
arable  to  the  thicknesses) .  The  region  where 
at**  f  becomes  comparable  with  unity  is 

x  -  a  =  0(k*)  =  0(e*) , 


indicating  that  the  leading-edge  singularity  wilt 
not  affect  the  usefulness  of  the  theory  until  one 
reaches  the  nose  region  which  is  x-a  »  0  (  t"*), 
since  -r «  0(«t  )  here. 

It  may  be  pointed  out  that  the  magnitude  of 
the  drag  coefficient  contributed  by  the  surface 
pressure  in  2-D  is  of  the  order  (of  which  a 

nonvanishing  portion  is  contributed  by  the  singular 
pressure  coefficient  near  the  nose  region),  while 
the  stagnation  region  at  the  nose  contributes  tc  a 
smaller  leading-edge  force  comparable  to  e*  ,  i.e. 

4*  .  The  latter  is  exactly  an  order  e  higher 
than  and,  therefore,  is  essential  if  the  3*D 

effect  on  the  transonic  drag  rise  is  to  be  consis¬ 
tently  analyzed,  for  this  reason,a  development  of 
the  lifting- line  theory  on  the  basis  of  the  octen¬ 
nial  eqimtiohw'i  thout  the  smal  I -disturbance  assumption 
appear  to  be  worthwhi le(5b) .  Additional  reasons 
favoring  the  use  of  a  full  potential  equation  are 
given  later  in  §  5. 


a  .portion  of_surface  containing  the  maxi  mum- veto- _ 

city  point,  with  the  exception  of  certain  very  spee- 
ial  distribution  of  f  (7) ■  Cook's  implicit .assumption  _ 
then  would  amount  to  the  stipulation  that  f (7)«?  ; 


Nonun i formi tv  Near  a  Shoe*  Root 
For  flow  with  an  imoeooed  supercritical 
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.'eg i on  terminating  by  a  shock,  the  reexpanding  sub* 
sonic  flow  develops  a  weak  singularity  at  the  shock 
root  (the  intersection  of  the  shock  and  an  imper¬ 
meable  surface),  unless  the  surface  has  a  zero  cur¬ 
vature,  as  is  well  knownlS".  In  particular,  the 
fluid  tangential  velocity  behaves  as 

(u-i ~  C*I.3«) 

/\  A 

where  u.  a  ft*  ,  the  subscript  "  sr  "  refers  to  the 
shock  root  and  £3  (+»*?;{  and  vp  are  rescaled 
coordinates  with  the  origin  a£  the  shock  root: 

7*  Zsrl&lJ'  £ 

A  %A 

i  n  above  Z*  «  is  proportional  to  the  sur¬ 

face  curvature  and  «c,  is  a  constant  determined  by 
the  surrounding  flow,  ft  follows  that  uuX  - 
becomes  infinity  at  the  shock  root  (f*o)  1  i ke 
leading  to  a  logarithmically  infinite  shock  jump 
for  £  3 $, /»x  at  the  shock  root,  (ef.  Eq.  (3.12a)). 
Although,  one  may  argue  that  this  rather  weak  sin¬ 
gularity  would  affect  little  the. lift,  moment  as 
well  as  drag,  the  failure  to  yield  a  correct  (fin¬ 
ite)  pressure  jump  at  the  surface  appears  to  have 
defeated  one  of  the  original  purposes  for  analysing 
the  shock  perturbation. 

However,  this  local  breakdown  can  be  treated 
(I kb),  once  if  it  is  realized  that  the_form  of  the 
reexpansion  singularity  in  u  (not  in  0)  given  in 
Eq.  (4-3)  is  unaffected  by  the  presence  of  the  3*0 
compressibility  correction  a®c  $a  y  in  the  composite 
POE  ( 3 - *♦)  -  Thus  any  3-0  (or  the  analogous  unsteady) 
influence  on  the  local  structure  of  the  if,  or  pres¬ 
sure  field, can  result  only  through  changes  in  the 
three  constants  ,  £sr  ,  and  .  In  fact, 

ay  expanding  uV|  *,r  and  k.  in  small  C  for  fixed 
x  and  z  (appropriate  at  points  far  from  the  shock 
root),  one  recovers  from  Eq.  (4.3a)  the  logarithmic 
singularity  in  the  previous  expansion.  Furthermore, 
latching  this  with  the  singular  solution  to  p  m 
$,  *  «  $,  +-  perrai ts  the  determination  of  the  correc¬ 
ted  values  for  the  shock  jump  and  for  the  shock- roct 
location. 

Implicit  in  the  formulation  involving  shocks 
-r.  |  3  is  that  c  3  $,  g  S  is  small  relative  to 
r&pD  •  This  condition  is  not  realized  where 
the  shock  strength  approaches  zero.  At  least  for 
solutions  constructed  from  the  similarity  structure 
'  §  3  -  7) ,  in  which  the  shock  jump  is  controlled  ex¬ 
plicitly  by  £d?/dy  ,  (weak)  unrealistic  expansion 
snock  may  apoear  from  such  an  analysis  (if  the  2-0 
solution  ‘or  the  component  flow  has  a  weak  shock). 
However,  the  inner  solutions  based  on  the  composite.. 
POE  (3-4),  or  on  POE  (3.26)  for  the  unsteady  an¬ 
alogy.  will  not  encounter  this  difficulty.  Compar¬ 
ison  pf  the  similarity  solution  with  that  based  on 
tne  unsteady  analogy  for  the  case  involving  weak 
snocxs  is  therefore  valuable. 

w. 4  Sonuni formi ties  at  Wing  Tips  and  at  the  Apex 


(4.3b) 


from  the  tips,  and  at  a  distance 

7  =  0(«) 


(4.5) 


from  the  apex  of  q  swept  wing.  Implicit  in  Eq. 
(4.4)  is  the  assumption  that  the  span  load  computed 
from  3o  (i.e.  from  the  strip  theory)  vanishes  at 
the  tip  like  an  elliptic  load  (ef.  Fig.  6a)>.  For 
an  untwisted  tapered  wing,  or  one  with  a  square  tip 
(cf.  Fig.s  6b  and  6c),  we  infer  that,  instead  of 
Eq.  (4.4),  the  nonuniformity  in  this  case  occurs  at 


y*i=0(«).  (u. 6) 

This  follows  from  the  last  integral  of  Eo.  (2.6c), 
which,  when  interpreted  in  Stieltjes1  sense,  gives 
the  additional  terms  (for  7-*0) 


* 


l  +  y 


[ I  *  SinA(-i)] 


(4.7) 


A  comparison  of  the  present  solutions  with 
more  exact  numerical  analyses  for  the  cases  in¬ 
volving  untwisted  tapered  swept  (or  oblique)  wings 
and  for  the  swept  wings  with  an  apex,  should  be 
helpful  in  assess inq  the  extent  to  which  these  non¬ 
uniformities  affect  the  usefulness  of  the  analysis. 

It  may  be  pointed  out  that  the  flow  field 
next  to  the  apex  must,  in  any  practical  situation, 
be  modified  by  the  presence  of  a  (slender)  fuse¬ 
lage;  therefore,  the  solution  breakdown  near  y  •  0 
Should  not  be  a  serious  practical  concern,  and  may 
presumably  be  separately  treated.  The  relatively 
small  scale  in  7  or  (y  fcl)  indicated  by  Eq.  (4.4)  - 
(4.6)  suggests  that  the  equation  governing  the 
regions  around  the  apex  and  the  tips  to  be  re¬ 
analyzed  is  the  3“D  TSC  equation  (noting  that 
£  3  J,'1’*  )•  These  extensions  remain  to 

be  investigated. 


5.  COMPUTATIONS:  EXAMPLES  AND  C0MPARISI0NS 
WITH  FULL  POTENTIAL  SOLUTIONS 

Inasmuch  as  the  existence  and  uniqueness  of 
the^solut ions  .cannot  be  easily  investigated,  dem- 
onstration  of  numerical solut ions_tp  the. reduced, 
problem  Is  an  essential  part  of  the  study.  The 
theory  has  been  iimited  by  the  assumptions  of  a 
high-aspect-ratio  and  of  the  small  disturbance;  it 
is  uncertain  that  the  analysis  may  predict  the 
aerodynamic  characteristics  for  swept  wings  to  the 
same  degree  enjoyed  by  the  lifting-line  theory  for 
the  straight  unyawed  wings.  In  view  of  this,  and 
the  several  nonun i formi t ies  of  the  analysis  dis¬ 
cussed  in  the  preceding  section,  we  consider  the 
direction  comparison  with  the  more  exact.  3*0  full 
potential  solution  an  integral  part  of  this  research 
Before  going  on  to  the  discussion  of  the  computa¬ 
tions,  some  remarks  will  be  made  on  the  basis  of  the 
more  exact  potential-flow  computation  programs  which 
are  involved  in  tne  comparison  studv. 


As  made  apparent  in  /2.3.  the  analysis  given 
in  §§  2  and  3.  without  further  refinement,  will 
preaxdown  at  a  distance 

/  i  I  -  O  ( 4.4) 


5. 1  Remarks  on  Existing  1-D  Potential  Computer 

Cod* 

The  unrefined  nature  of  the  far-fleld  descrip¬ 
tion  in  most  discretized  flow-field  computation 
methods,  as  a  result  of  the  computer  storage  and 


'0 


computation  time  limitations,  >s  well  Known,  "his 
problem  becomes  more  serious  in  the  3*0  cases  .  _ 
(20*24).  It  is  not  at  all  clear  from  the  published 
data  whether  the  grid  distributions  used  therein 
are  sufficiently  refined  for  the  purpose  of  de* 
scribing  adequately  the  upwash  induced  by  the  far- 
wake  vortici ties, which  has  been  so  crucial  to  the 
analysis  of  a  high-aspect-ratio  wing. 

There  appears  to  be  an  additional  problem 
brought  about  by  the  scarcity  of  the  span  stations 
available  in  the  current  3*0  transonic  program 
(21-24),  jn  tne  FLO  22  code  for  planer  wings '22) 
a  total  of  twenty-one  (21)  span  stations  are  al¬ 
lowed.  Applying  it  to  an  oblique  wing,  for  example, 
there  remains  only  ten .(10)  stations  on  each  wing 
pane!,  and  one  certainly  cannot  attach  too  much 
confidence  to  the  results  obtained  for  the  tip 
region  in  this  case.  It  has  also  been  known  among 
users  that  the  span  load  and  sectional  lift  co¬ 
efficient  so  obtained  are  far  from  being  very 
smooth— the  values  at  successive  stations  appears 
to  alternate  noticeably  about  some  mean  curve, 
fig.  7  shows  typical  results  of  sectional  lift  co¬ 
efficient  obtained  from  the  converged  solutions 
generated  from  two  versions  of  the  FLO  22  code  for 
ah  'elliptic  wing  pivoted  at 30?at  H  —  »  0.755.  Except 
for  this  limited,  but  noticeable,  irregularity  in 
the  transonic- speed  range  which  require  some 
caution  in  making  a  comparison  with  the  asymptotic 
analysis  the  method  is  known  to  be  capable  of 
reproducing  the  span  loads  determined  by  the 
(linear)  pane!  methods  for  the  subsonic  speed 
range (58) . 

Presumably,  this  irregularity  has  not  been 
observed  in  the  3-0  computations  using  the  tran¬ 
sonic  smai 1 -disturbance  equation'20).  We  may  point 
out  however  that  the  current  3-D  TS0  codes  are  not_ 
applicable  direct) Ytoproblems  lacking  *  bi-lateral 
symmetry  such  as  that  of  an  oblique  wing.  _Futher-_ 
.more,.  .  _the  smal '-disturbance  assumption  breaks 
down  near  the  leading  edge,  comparison  with  TSD 
code  will  not  reveal  one  of  the  imoortant  limita¬ 
tions  of  the  present  ( transonic )  theory.  For  the 
two  above  reasons,  comparison  with  the  TSD  codes 
has  not  been  made.  We  believe  such  a  comparison 
can  be  quite  useful  in  our  future  study  focusing  on 
shock-wave  patterns. 

From  a  strict  theoretical  view  point,  the  most 
serious  limitations  of  current  3-0  full-potential 
codes  is  the  empiricism  introduced  by  modelling  the 
inviscid  wake  in  the  formulation,  which  we  believe 
is  unnecessary.  Two  assumptions  were  introduced 
in  Refs.  21-24:  (!)  the  shape  of  the  trailing 

vortex  (TV)  sheet  is  speci fled  a  priori ;  (ii)  a 
condition  equivalent  to  (in  the- present  notations) 

<PMX=  O  (5.1) 

is  applied  at  the  TV  sheet.  Assuming  a  thin  air¬ 
foil  section,  assumption  (i)  is  seen  to  be  equiva¬ 
lent  to  that  in  the  TSD  theory,  with  the  TV  sheet 
transferred  to  the  "wing  trace"  on  the  x-exls.  The 

condition  (ii)  is  incorrect,  strictly  speekinq,, _ 

and  requires  a  clarl f ication  .  Were  Eq.'TS.  IT 
to  hold  on  both  sides  of  the  sheet,  Refs.  22-24 
would  imply  that  the  FOE  at  far-downstream  should 
yield  (in  the  present  notations) 


*  This  controversial  condition  was  brought  to  my 
attention  by  Or.  Norman  0.  Malmuth. 


o  .  or  5.31 

that  is,  the  upwash  is  continuous  across  tne  sneet 
or  the  trace.  This  and  the  continuity  of 
are  then  the  correct  TV-sneet  condition,  wnicn  is  at 
least  consistent  with  the  smal I -disturbance  theory. 

On  the  other  hand,  the  programs  similar  to  those  in 
Refs.  22-24  do  lead  to  results  consistent  witn  ex¬ 
perimental  data.  This  can  be  reconciled  by  tne 
observation  that  if  we  were  to  stipulate  the  in¬ 
correct  assumption  of  continuity  in  ,  then  Eg. 
(5.3)  would  become  consistent  with  ^  *  0.  it  is 
apparent  that  in  the  difference  equations  of  Ref. 
22-24,  Eq.  (5.3),  but  not  Eq.  (5-1).  was  actually 
used.  We  may  conclude  that  whilej^w  0  appears  to 
be  erroneous,  the  code  actually  uses  the  more 
appropriate  condition  Eq.  (5-3),  and  that  the  wake¬ 
modelling  in  Ref.  21-24  introduces  an  error  com¬ 
parable  to  that  in  the  smal I -disturbance  theory. 

5.2  Computations  of  9.,  and  R 

The  reduced  mixed-type  problems  for  the  , 

6,  ,  and  in  the  similarity  solutions  (  §  3-7) 

are  solved  numerically  by  a  relaxation  method,  using 
type-sensitive  difference  operators  corresoording  to 
Wurman' s  "fully  conservative  form’';P).  The  pro¬ 
cedure  of  line  relaxation  used  for  may  be 
considered  standard,  except  for  the  uses  of  an  Im¬ 
proved  far-field  description  with  Eq.  (3.23a)+  and 
of  a  third-order  convergence  acceleration  scheme/60) 
Shock  fitting  aigori thm'6l )  js  being  used  in  con¬ 
junction  with  the  sjnock  perturbation  analysis  based 
on  Eq.  (3.22)  for  0,  ,  and  similar  treatment  for  , 
but  will  not  be  discussed  here,  inasmuch  as  the 
cases  with  subcritical,  and  slightly  supercritical, 
component  flow  studied  below  in  S  5-3  are  shock- 
free. 

For  ,  and  $4  ,  the  problems  are  linear 
and  the  transition  and/or  shock-jump  boundaries  are 
known  from  the  solution,  but  they  require  a 
greater  storage  owing  to  the  variable  coefficients. 

A  nonuniform  mesh  with  81  x  61  grid  points  over  a 
computation  domain  It 'I  £  6,  l£l  £  6,  were  used  bv 
Cheng  and  Meng('3,l4)(  is  uniform  over  the  wing 

section).  A  subsequent  program  employs  three 
successive  grid-halvings  up  to  a  mesh  of  124  x  32 
(with  nonuni  form  aZ  on  the  wing).  The  leading  and 
trailing  edges  are  located  at  x  ■  -’  and  x"  »  +1, 
respectively  (l.e.  -  a  *  b  »  I);  a  departure  of 
from  negative  one  can  be  accounted  for  by  changing 
all  x  to  (x- ♦  I  +  a)  in  the  equations  governing  the 
similarity  solutions.  The  case  -a"  ■  •  I  cor¬ 

responds  to  a  swept,  or  oblique,  wing  with  the 
straight  axis  x  »  0  coinciding  with  the  50*.  chord 
i  i  he . 

The  iterative  solutions  in  Refs.  13  and  14  uses 
uses  a  relaxation  factor  of  1.8  and  0.8  in  the 
elliptic  and  hyperbolic  region,  respectively;  200- 
300  sweeps  are  needed  for  convergence  of  the  circu¬ 
lation  to  within  10*5.  On  an  IBM  3031,  using 
double-precision  arithmetics,  the  convergenced 
solution  for  each  ^  requires  typically  10-12 
minutes*  We  note  in  passing  that  a  computer 

*  A  2-5%  change  in  surface  pressure  results  from 
the  Improvement  in  the  far-field  description 
based  on  Eq.  (3.23a). 

*  The  IBM  3031  at  tne  University  of  Southern 
California  Campus  has  an  expandable  core  memory 
of  300-6400  K  words . 


program  similar  to  that  for  hat  baan  cons i da rad 
in  the  context  of  a  straight  unyewed  wing  by 
Smell*1''. 

5. 3  Examples:  Comparison  with  Solutions  Band 

on  tha  Ful  l-Fotandal  ROE 

for  tha  expressed  purpose  of  comparing  tha 
asymptotic  analysis  (tha  similarity-solution 
structura  in  particular)  with  tha  corrasponding 
ful l-potantial  solutions,  wa  considar  below  exam- 
plas  of  obliqua  as  wall  as  symmetric  swept  wings 
with  high  subcritical  and  slightly  suparcri tlcal 
conponant  flows.  All  wing  pi  an forms  consldarad 
ara  alllptic,  with  tha  major  axis  coinciding  with 
tha  mid  (502)  chord;  tha  wing  s act  ions  ara  ganar- 
atad  from  a  single  profile  NASA  3612-02,  60, 
rescalad  to  an  arbitrary  thickness  Tha  latter 
will  be  sat  equal  to  d  in  tha  theory.  Thesa  wings 
fulfill  tha  geometrical  description  of  £q.  (3.18) 
and, therefore,  tha  inner  solution  can  be  obtained 
by  a  I i near  combination  of  tha  basic  similarity 
solutions  9i  ,  and  &  ,  which  will  be  de¬ 

termined  after  specifying  tha  component  transonic 
parameters  K.  „ ,  the  i nc i dance , and  tha  locations  of 
tha  leading  and  trailing  adgas.  Via  pointed  out 
that  the  airfoil  profile  considered  has  bean  used 
in  various  wind  tunnels  and  preliminary  design 
studies  of  obliqua  wings  at  H  ■  0.60  -  1.4(6,62). 
Several  sets  of  these  basic  similarity  solutions 
have  been  obtained  fqr  K  „  •  3.6  and  Kn  “3-65, 
and  have  been  described  in  tone  data  1 1  in  Refs .  13 
and  14.  These  solutions  are  used  in  the  subse¬ 
quent  comparison  with  the  full  potential  solutions. 

To  illustrate  tha  solutions*  behavior,  we 
jjjow  tha  surface  distributions  of  i0£,  T,g,  and 
a,-  in  Fig.  3.  for  zero  incidence  jnd  Kn  •  3-65. 

The  basic  component  flow,  for  the  i0jr  shown  in  Fig. 
3  is  siightly  supercritical.  The  critical  condi¬ 
tion  for  the  2-0  componen t  in  this  casa  Is 
fox  *  *"/(*♦  1)  -  I.43&. 

Numerical  results  comparable  with  our  soju-  _ 
tions  are  generated  from  one  version  of  A. 

Jameson's  3-0  full-potential  computer  codes 
"FLO  22"  (See  Refs.  6,21),  which  is  used  with  im¬ 
plementations  for  oblique-wing  analyses  at  NASA 
Ames  Research  Center  Aeronautics  Division  and  at 
Grumman  Aerospace  Corporation  Research  Department. 
The  algorithms  employed  in  FLO  22  are  not  fully 
conservative,  but  this  may  not  ba  essantial  for 
shock-free  solutions  presented  below.  We  point 
out  that  the  rLQ  22  data  from  NASA  Ames  and  from 
Sruiman  are  not  identical,  owing  mainly  to  tha  use 
of  different  meshes.  The  availability  of  data  . 
from  two  sources  is  helpful  in  delineating  the  _ 
nature  of  qiscreoancy  between  our  theory  and  cue 
more  exact  3-0  programs  as  noted  in  §5.1.  Oata 
from  the  latter  is  still  influenced  by  the  mesh 
site,  spacing  of  the  span  stations,  number  of  It¬ 
erations,  the  detail  of  the  leading-edge  geometry 
description,  which  are  different  in  the  NASA  Ames 
and  the  Grumman  calculations. 

A  number  of  FLO  22  runs  have  been  made  for  the 
obiiaue  wings  with  free-stream  Hach  number,  swept 
angle,  wing-thickness,  etc.  chosen  to  give  either 
<n  •  3.6  or  Kfl  «  3.65  employing  the  same  basic 
airfoil  section.  An  elliot’c  planform  is  used  in 
saen  case;  wing  twist  and  wing  bend  are  assumed  to 
pc  zero,  as  noted  earlier.  Among  the  first  com¬ 
parison  studies  made  is  a  case  with  relatively 
thin  wing  section  (64  thickness  ratio)  and  a  rath¬ 
er  nign  ascect  ratio  (an  elliptic  planform  with  a 


major-to-minor  axes  ratio  of  20).  The  surface  pres¬ 
sure  coefficients  in  this  case  nave  been  presented 
in  the  journal  paper  (Ref.  !4e),  we  consolidated 
theei  here  in  a  single  plot  (Fig.  5)  in  order  to 
show  differences  from  the  subsequent  comparisons 
involving  considerably  thicker  wing  sections  and 
larger  sweep  angles.  The  free-stream  Mach  number 
and  the  sweep  angle  used  In  the  FLO  22  calculation 
is  0.8242  and  A»  22.5  ,  giving  a  component 
Hach  number  H  ■  0.7615.  Thus,  we  have  Kn  ■  3.60, 
0-  1.003  and  <  •  0.1277  in  this  case;  the  com¬ 
ponent  flow  is  slightly  below  being  critical.  The 
FLO  22  data  from  NASA  Ames  (in  ttssl  I  crosses  and 
"v*'J  and  from  Gr unman  (In  smell  open  circles)  shown 
in  the  plot  appear  to  be  rather  close  except  next 
to  the  leading  edge.  The  Cp  values  computed  from 
the  similarity  solution  (shown  in  solid  curves) 
agree  reasonably  well  with  the  FLO  22  data  except 
near  the  leading  edge,  as  anticipated. 

Encouraging  is  that  the  degree  of  agreement 
with  the  FLO  22  data  does  not  appear  to  deterio¬ 
rate  much  with  increasing  the  wing  thickness,  or 
reducing  the  wing  aspect  ratio.  The  consolidated 
plot  in  Fig.  10  shows  the  surface  pressure  coeffl  - 
dent,  at  seven  span  stations  on  a  122  thick,  14:1 
elliptic  wing,  pivoted  at  A  *  30  for  H— ■  0.7677. 
The  coaiponent  Hach  number  is  Hn  -  0.6648  in  this 
case,  giving  K„  -  3.45.  The  three  sets  of  surface 
dete  of  Fig.  8  are  used  to  construct  the  Cp  distri¬ 
butions  shown  in  heevy  solid  curves.  The  Cp  value 
for  the  componen t  flow  is  -0.689  in  this  case  and 
both  solutions  give  the  appearance  of  supercritical 
shock-free  regions ton  the  upper  surface  of  the  down¬ 
stream  wing  panel.  ' 

The  agreement  between  the  FLO  22  data  and  those 
based  on  the  asymptotic  analysis  in  both  Figs.  9 
and  10,  should  be  considered  as  being  better  than 
expected,  inasmuch  as  the  relative  error  in  the 
asymptotic  theory  belongs  to  an  order  determined  by 
or  C*.  which  ever  Is  larger.  The  magnitudes 
of  for  the  examples  shown  in  Fig.  9  and  Fig. 

10,  are,  respectively ,  *ev*  ■  0.153  and  0.243.  it 
may  be  recalled  that  there  Is  a  noticeable  differ¬ 
ence  between  the  two  sets  of  FLO  22  data  shown  in 
Fig.  7.  The  eg remnant  found  here  appears  to  add 
credence  to  the  Ames  data  set  in  this  case. 


I  t~is  not  altogether  clear  that  the  ancourag 
ing  agreement  found  for  the  pivoted  wings  may 
still  hoid  to  some  degree  with  the  presence  of  an 
apex  in  a  symmetric  swept  wing.  Figure  II  presents 
a  consolidated  plot  for  surface  pressure  coeffi¬ 
cient  at  seven  span  stations  on  a  synsnetric  swept 
wing  with  the  same  basic  (elliptic)  pi  an form, sweep 
angle, and  section  profile  as  those  in  the  preceding 
ing  figure.  The  free-stream  Hach  number  is  however, 
lower  with  H—  •  0.7549;  thus,  H„  .  0.6538, 

■  J;6ff;  •  ■  9.  tdd*  an«  ®  •  I . W2.  The  sot ut Ton 
remains  subcritical  for  most  stations,  component- 


"In  both  the  asymototic  analysis  and  the  FLO  22  cal- 

dicular  to  the  straight  axis,  therefore,  the  sec¬ 
tional  cuts  shown  in  Fig.  to  are  graphically  incor¬ 
rect.  The  differences  resulting  from  changing  the 
orientation  of  the  cuts  turns  out  to  ba  numerically 
small,  however. 


co  illustrate  tut  manner  in  wnicn  tne  solutions  be¬ 
have  nur  the  tip,  wnera  ch«  t rands  of  departure  of 
(ha  asymptotic  from  tha  FLO  12  ana  Witt ,  at  well  at 
between  tha  two  tats  of  FLO  22  data  themselves,  art 
claariy  avidant.  Away  from  tna  tip,  agraamen t  it 
at  rpasonaala  at  in  tha  two  preceding  caaiparitont 
for  tha  oblique  wings.  However,  tha  tat  of  FLO  22 
data  in  fill m  c ire  tat  (black  dots)  g i vat  a  conti r 
tartly  highar  paak  for  -Cp  than  tha  ona  in  opan  cir- 
clat.  Tha  ditcrapaney  aay  da  explained  dy  tha  un¬ 
carts  i  nty  atsociatad  with  tha  fluctuating  tpanwiaa 
dittridut ions  of  tactional  lift  shewn  In  Fig.  J. 

Tha  asymptotic  analysis  lands  to  giva  a  conslttantly 
I owar  -Cp  than  tha  FLO  22  data  on  tha  lower  surface 
around  tha  quarter  chord.  This  tendency  is  also 
apparent  froa  Fig.  10.  This  saw  1 1  but  noticeable 
discrepancy  way  be  associated  with  tha  accuracy  and 
degree  of  convergence  in  the  11  solution.  Results 
obtained  moat  recently  for  Tq  based  on  a  wort  re* 
fined  wash  and  with  a  saw  1 1  residua  appears  to 
alter  tna  Cp  values  on  the  lower  surface  slightly 
in  the  right  direction. 

The  nost  important  piece  of  inforawtion  f row 
Fig  II  it  the  coaipar i son  made  for  the  station 
closest  to  the  apex.  At  I0R  semi -span  from  tha 
tyaewtry  plane  (y  •  0.097).  the  agreement  of  tha 
asymptot ic  results  with  the  FLO  22  data  (In  opan 
circles)  raawin  as  good  as  other  stations.  Evan 
for  the  FLO  22  data  belonging  to  the  second  tat 
(black  dots) ,  the  agreement  with  tha  results  it 
on  the  similarity  solution  is  still  reasonable. 

6.  CONCLUDING  REMARKS 

Tha  foregoing  presentation  has  shown  that 
Frandtl's  lifting-line  idea,  originally  applied  to 
an  unyawed,  straight  wing  of  high  aspect  ratio,  can**" 
be  extended  to  tha  study  of  3-0  mixed  flows  over 
transonic  swept  wings.  Examples  of  oblique 
(pivoted)  and  symmetric  swept  wings  involving  high 
tubcritical  and  supercr I t i ca I  component  flaws  ere 
shown;  comparisons  with  full-potential  solutions 
from  existing  codas  are  made.  Except  near  the  wing 
tip  and  tha  wing  root,  and  also  tha  leading  edge, 
where  breakdowns  are  expected;  the  agreement  with 
-the  full-potential  results  i a  encouraging.  j 

A  salient  feature  brought  out  by  the  asymptot¬ 
ic  analysis  is  a  contribution  to  the  upwath  correc¬ 
tion  proportional  to  ft mtfQtf1/*? ,  resulting  from 
the  spanwise  (y1-)  component  of  the  locally  chadded 
wake  vorticities,  even  though  the  latter's  affect 
on  tha  other  velocity  component  is  relatively 
small  (in  the  transonic  regiam  considered).  The 
dependence  of  the  upwath  correction  on  the  sweep 
and  on  the  span  loading  have  been  shown;  the  expli¬ 
cit  result  presented  in  f  2.1  should  be  quite  use¬ 
ful  for  controlling  the  3-0  effects  on  the  transon¬ 
ic  component  flow  in  design  studies.  Another  at  -  \ 
sent i el  ingredient  of  the  analysis,  absent  in  the 
classical  theory,  is  a  source-ilka  coa^ressibi lity 
correction  to  tha  FOE  of  tha  component  flow;  it 
leads  to  a  line-source  effect  in  tha  outer  flow. 

An  important  gain  in  tha  theoretical  develop¬ 
ment  through  this  approach  is  the  availability  of 
a  similitude  in  the  3*0  flow  structure  (  fi.7), 
applicable  to  wings  with  sections  generated  from  a 
single  profile.  With  the  similarity  solutions, 
the  reduced  2-0  problems  can  be  solved  once  for 
all  span  stations;  the  explicit  dependence  of  the 
flow  field  on  wing  sweep,  aspect  ratio,  and  span- 
wise  distribution  of  the  local  chord,  as  well  as 
twist  and  band,  can  be  studied. 


Tha  examples  analyzed  tnus  far  do  not  include 
cases  with  imoedded  shocks.  Solutions  making  use 
of  tha  3-0  corrections  for  the  shock  jump  and  the 
shock  geometry  (  $ 1 3 . 3  and  3-7)  and  of  the  treat¬ 
ment  of  the  re-expension  nonuniformi  ty  I  J  k.  3) , 
remain  to  be  demonstrated.  In  this  paper,  we  nave 
not  discussed  the  computational  study  based  on  tna 
unsteady  analogy  <  f  3-&) ,  yet  to  be  completed.  Un¬ 
published  results  computed  from  this  analvogy  with 
an  ADI  algorithm  by  Evans  compare  rather  closely 
with  similarity  solutions  and  promise  an  effective 
procedure  for  treating  complicated  shock  patterns. 

Whereas  a  localized  treatment  of  the  leading- 
edge  breakdown  is  possible,  a  solution  to  this 
difficulty  is  to  jtevelop  a  lifting-line  theory  on 
tha  basis  of  the  potent  ITT  theory  without  the  small- 
disturbance  assumption.  A  work  along  tnis  line  net 
already  been  carried  out  in  Ref.  58,  where,  inter¬ 
estingly,  solutions  via  similarity  flow  structure 
may  again  be  found.  As  noted  earlier,  in  (12.1), 
the  inner  solution  in  such  a  formulation  should 
provide  a  proper  treatment  of  the  vortex-sheet 
geometry  near  the  trailing  edge,  unaccounted  for  in 
most  (if  not  all)  current  3-D  computer  codes. 
Treatamnt  of  nonuniformities  at  the  root  and  the 
tip  sides  are  the  remaining  problems  for  the  lift¬ 
ing-line  theory;  their  analyses  for  the  transonic 
regime,  as  observed  in  (  Jk.k)  .would  involve  solv¬ 
ing  a  fuller  3~D  problem. 
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FIGURE  CAPTIONS 

Fig.  I  The  bound  and  trailing  vortex  system  il'us 
strated  for  (a)  a  pivoted  (ooiique) 
straight  wing,  lb)  swept  wingfwith 
straight  center-line  segments,  and  (c) 
swept  wings  with  curved  center  line._ 
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Fig.  2  The  Cartesian  end  the  orthogonal  eurvlllin- 
Mr  coordinates  illustrated  for  the  wing 
plan*.  Tha  fc  and  * 'ordinates  ara  normal 
to  tha  wing  plana. 

Fig.  J  Illustration  of  a  spanwisa  component  of  cha 
trailing  vorticity  which  indues*  to  a  log* 
archaically  infinita  upwash. 

Fig.  4  Tljn  upwash  function  I II us t rated  for 

8.39  at  three  yaw  angles  for  an  el  I  ipn 
tic  end  an  extended-span  distribution  In  f* 


Fig.  S  Tha  upwash  function  £*  2£l  I  lustrated  for 
el  Optically  loaded  swapt-forward  and 
swept-back  wings  of  ^  ■  14  at  N—  » 
0.755.  Refer  to  text  for  conversion  to 
other  aspect  ratio  and  Mach  numbers. 

Fig.  6  Sketches  Illustrating  different  types  of 
wing  tip  geometry:  (a)  elliptic  or  para¬ 
bolic  type,  (b)  tapered  swept-wing 
(standard)  (c)  tapered  swept  wing. 

Fig.  7  Spanwisa  Distribution  of  sectional  lift 

coefficient  computed  for  an  oblique  wing 
by  solutions  to  the  full  potential  equa¬ 
tion  (FLO  22)  from  two  sources. 

Fig.  8  Similarity  solutions  for  * 'of-  fig.  and  f2J? 
on  the  upper -and  lower  surfaces  of  an 
oblique  wing  at  Kn  *  3.45.  The  straight 
axis  is  located  at  the  mid  chord;  the  air¬ 
foil  section  Is  generated  from  NASA  3812- 
02,40  rescaled  to  an  arbitrary  thickness; 
the  section  Is  set  at  zero  Incidence,  with¬ 
out  twist  and  without  wing  bend. 

Fig.  9  A  consol idatad  plot  showing  surface  Cp  on 
three 'span' stations  of  a  hTgh-espect-  ra¬ 
tio  wing  pivoted  at  22.5*  ,  at  JFJ  ight 
Mach  number  9. 82k.  The  planform  is  a  29:1 
ellipse;  the  wing  section  is  NASA  3806- 
02,40  set  at  zero  incidence,  without 
twist  and  without  bend. 

Fig.  10  A  consolidated  plot  showing  surface  Cp  on 
seven  span  stations  of  an  oblique  wing 
pivoted  at  30*  ,  at  Flight  Mach  number 
0.767.  The  planform  is  a  54:1  ellipse; 
the  wing  section  is  NASA  3612-02,40,  zero 
incidence,  zero  twist  and  zero  wind  bend. 

Fig.  II  A  consolidated  plot  snowing  surface  Cp  on 
seven  span  stations  of  a  symmetric  swept- 
aeck  wing  with  30*  sweep  angle,  at  Flight 
Mach  number  of  0.755.  The  basic  planform 
is  a  14:1  ellipse;  the  wing  section  is 
NASA  3612-02,40,  zero  Incidence,  zero  twist 
and  no  wing  bend. 
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